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Abstract 



We calculate saturation curves for strange hadron matter using recently con- 
structed baryon-baryon potentials which are constrained by SU(3) symmetry. 
All possible interaction channels within the baryon octet (consisting of N, 
A, S, and H) are considered. It is found that a small A fraction in nuclear 
matter slightly increases binding, but that larger fractions (> 10%) rapidly 
cause a decrease. Charge-neutral {N, A, H} systems, with equal densities for 
nucleons and cascades, are only very weakly bound. The dependence of the 
binding energies on the strangeness per baryon, f s , is predicted for various 
{N, A, H} and {N, A, systems. The implications of our results in rela- 

tivistic heavy-ion collisions and the core of a dense star are discussed. We 
also discuss the differences between our results and previous hadron matter 
calculations. 
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I. INTRODUCTION 



The study of the properties of strangeness-rich systems is of fundamental importance 
in understanding relativistic heavy-ion collisions (1J and some astrophysical problems 0. 
The qualitative features of such systems and their possible detection in the universe and in 
relativistic heavy- ion collisions were first discussed by Bodmer [|3] in 1971. Within Quan- 
tum Chromodynamics (QCD), it was suggested that the strangeness-rich systems could be 
strange quark systems consisting of up (w), down (d), and strange (s) quarks. These exotic 
systems could be either metastable states || against the decays into hadrons, or absolute 
bound states with energies much lower than normal nuclear matter |]||. However, the theo- 
retical calculations for strange quark systems are still in the developing stage. For example, 
within the MIT bag model it is found J7j] that the stability of the bag strongly depends 
on the rather uncertain bag constant -Bbag- The strange quark matter is absolutely stable 
for ~ 140 MeV, metastable for = 150-200 MeV, and unstable for > 200 
MeV. Therefore, one cannot rule out the possibility that the strangeness-rich systems could 
be strange hadron systems made of nucleons and hyperons; as studied, for example, in 
Refs. gHTJ. 

The production and detection of strangeness-rich matter from relativistic heavy-ion col- 
lisions has been studied in recent years JT5^[T7]J7[] . Two scenarios have been discussed. The 
first one is the coalescent mechanism |17], which assumes that the produced hyperons are 
captured by the nearby nuclear fragments in the freeze-out region to form multi-hyperon 
clusters. The second mechanism is the distillation process [|l3|-|i5|,|7|j associated with the pro- 
duction of a quark-gluon plasma (QGP) in the baryon-rich region. The essential idea is that 
the s quarks in a QGP, in which ss pairs are abundant, are captured by the surrounding u 
and d quarks, liberated from the initial heavy ions, to form K + and K°. The emission of 
these kaons and other mesons causes cooling of the QGP into a strange quark system. The 
interesting question to ask is whether this system will be absolutely stable or metastable 
against the collapse into a hadron system. This question can be answered by comparing 
the energy of such a system with the energy of a strange hadron system with the same 
strangeness quantum number. It is therefore important to develop theoretical approaches 
to calculate the energies of strange hadron systems. In this paper, we report on the first 
results of our efforts in this direction. The calculation for the strange quark systems with 
similar sophistication is beyond the scope of this paper. 

Most of the recent investigations of strange hadron systems have been done by using the 
relativistic mean-field model P-[Tl"|,|7|1 . In addition to the usual a and u mesons, these models 
also contain a* and mesons, introduced in order to have strong attractive hyperon-hyperon 
interactions. The vector coupling constants are chosen according to SU(6) symmetry, while 
the scalar coupling constants are fixed to hypernuclear data. Extensive calculations for the 
systems consisting of {p, n, A, S°, H~} mixtures have been performed. Clearly, these calcula- 
tions are not completely consistent with the SU(3) symmetry, since S's are not included (for 
reasons which will be discussed below). Furthermore, it is not clear that the values for the 
coupling constants employed in these models are consistent with the very extensive data on 
nucleon-nucleon (AW) and hyperon-nucleon (YN) reactions. A rigorous prediction should 
be consistent with both the two-body data and the data of hypernuclei. 
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An alternative approach is based on the many-body theory with baryon-baryon potential 
models. This was first pursued by Pandharipande || using a variational method and rather 
crude baryon-baryon potentials. This approach has recently be revived in Ref. by using 
the Brueckner-Hartree-Fock approximation and the Nijmegen soft-core YN potential [18] of 
1989. There, the authors only consider an infinite system consisting of A's and nucleons. It 
was found that for a given total baryon density ps the binding energy per baryon, Eb = 
—E/Ab, decreases as the fraction of strangeness, f s = \S/Ab\, increases. The most stable 
system in their approach has Eb = —E/Ab ~ 16.0 MeV, at f s ~ 0.06 and p ~ 0.23 fm~ 3 . 
Comparing with the calculations using relativistic mean-field models as described above, this 
investigation is rather incomplete since the role of S is not explored, owing to the restriction 
of the employed Nijmegen potential. The importance of S was pointed out in Ref. ||. It 
is needed to stabilize the system against the strong A + A — > H + iV process, which can 
occur at relatively low density of A since the threshold energy for this reaction to occur is 
only about 28 MeV. The presence of H will Pauli block this reaction. Furthermore, their 
calculations do not include hyperon-hyperon interactions and hence the effects of additional 
hyperons on the hyperon mean field are neglected in solving the self-consistent G-matrix 
equation. In this paper we try to be as complete as possible in that we include all possible 
interaction channels that are allowed for the baryon octet. We use the recently constructed 
NSC97 baryon-baryon potential models |I^,2"0]j to describe all these channels. 

The content of the paper is as follows. In Sec. |H] we briefly highlight some of the features 
of the employed baryon-baryon potential models. In Sec. [I f| we review the definition of 
the G matrix and related quantities, with some emphasis on the treatment of the coupled 
channels. In Sec. |TV] we present and discuss the results of our calculations. We conclude 
with a brief summary of our findings in Sec. [V|. 



II. BARYON-BARYON POTENTIAL 

In this paper, we report on the first results from an investigation of strange hadron 



matter using the most recently developed baryon-baryon potentials [19 , 20| . These poten- 
tials are constructed within the dynamics defined by the SU(3) symmetry, and the data of 
nucleon-nucleon (AW) and hyperon-nucleon (YN) reactions. We follow a similar Brueckner- 
Hartree-Fock formulation as recently employed by Schulze et al. ||12|| , but with an important 
improvement: the starting coupled- channel AW © YN © YY potentials include the AA, AE, 
EE, and SA" channels of strangeness S = —2, as required by SU(3) symmetry. The presence 
of E then also enforces us to consider the SA and SE channels of strangeness S = —3, and 
the SS channel of strangeness S = —4. The assumption of SU(3) symmetry allows us to 
unambiguously define the baryon-baryon interactions for the S = —2, —3, —4 systems from 
the previously constructed AW and YN interactions. 

However, because of the lack of sufficiently accurate YN data and some uncertainties 
in SU(3) coupling constants, the constructed baryon-baryon potentials have some model 
dependence. In Ref. |[19| , six YN models have been constructed, based on different choices 
for the vector-magnetic F/(F + D) ratio, a™. Values range from a™ = 0.4447 for model 
NSC97a to a™ = 0.3647 for model NSC97f. The different choices for oiy (consistent with 
static or relativistic SU(6) predictions) are chosen such that the models encompass a range 
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of scattering lengths in the SiV and AA" channels, but all models describe the YN scattering 
data equally well. Differences show up in more elaborate applications such as hypernuclear 



calculations; see Ref. [B5 for a more detailed discussion. In this paper we will only con- 
sider models NSC97e and NSC97f, which seem to be the most consistent with the existing 
hypernuclear data JTTjj . 

Another important assumption in the construction of the baryon-baryon potentials is 
that the SU(3) symmetry is applied to the full range of the interaction; i.e., to the long- 
range as well as to the short-range part. Although there is no empirical evidence that the 
short-range part indeed satisfies the SU(3) symmetry (there are no YY scattering data to 
test this assumption, for example, with these potential models), we have chosen for this 
approach since it allows us to extend the AW and YN interactions to all YY interactions 
describing the S = —2, —3, —4 systems, without having to introduce any new parameters. 
To illustrate the differences between the various interactions, we show in Fig. [I] the 1 Sq 
elastic phase shifts for the two models in the NN, AA, EE(T = 2), and HS channels. The 
differences between NSC97e and NSC97f are fairly small and, at this scale, will only show 
up in the AA channel. A more detailed description of the NSC97 potential models in the 
YY channels will be presented elsewhere [pO . 



III. G-MATRIX FOR COUPLED CHANNELS 

In the next section we will present the results of several Brueckner-Hartree-Fock calcu- 
lations using the baryon-baryon potentials for the S = 0, . . . , —4 systems, and so here we 
review some of the aspects of the G matrix and define the relevant quantities. In each case, 
the (strange) nuclear matter is characterized by a total density p, which is broken up into 
the contributions from the four baryon species according to 

P = Pn + Pa + Ps + Ps 
= p(Xn + Xa + Xe + Xs) 

= JL ( 2 A^ 3 + 4 A)3 + 3k^ 3 + 2kf )3 ) . (1) 

This also defines the Fermi momentum k F for a baryon B with density fraction xb- As 
standard, we define the G matrix G(p[, p f 2 ; pi, P2) for incoming (unprimed) and outgoing 
(primed) momenta. Defining mass fractions pi = M^j \M\ + M 2 ), the total momentum is 
P = Pi + P2 = Pi + P2, and the relative momenta are given by k = P2P1 — P1P2 and 
k' = P2P1 — PiP2- The G matrix satisfies the Bethe-Goldstone equation |21 

G(k', k; P, u) = V(k', k) + [(%V(k', q) — p r G(q, k; P, w), 

J (2ny uj - £?i(piP + q) - E 2 (p 2 P-q) 

(2) 

where Q(q, P) is the Pauli operator which ensures that the intermediate- state momenta are 
above the Fermi sea (see below). The uj denotes the starting energy (including rest masses), 
while the intermediate-state energies are given by 
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1/2 112 

Ei(j® + E 2 {v>i) = M i + M 2 + 7^ + 7^f 2 +Ke U{p'{) + Tie U(j® 

= M+ ^M + Y^ +1Ze f/ ^ lP + q ) + Ue U(jMP-q), (3) 

with M and fi the total and reduced mass, respectively. The single-particle potentials U are 
defined by 

U(pi) = J j^G[p 1 ,p 2 ;p 1 ,p 2 ;u = E 1 (p 1 ) + E 2 (p 2 )}. (4) 

Hence, we have two equations, Eqs. (0) and (f|), which have to be solved self-consistently. 

In our calculations we have made several approximations. First of all, the energies 
are treated in the nonrelativistic expansion, as is obvious from the 1/M 2 expansion in 
Eq. The same 1/M 2 expansion was used in the derivation of the NSC97 potentials []T9 



Brockmann and Machleidt 53] have demonstrated the effect on nuclear-matter results when 
one uses instead the relativistic energies and the Dirac equation for the single-particle motion; 
these type of calculations have become known as Dirac-Brueckner calculations. For purely 
nuclear matter they find that the saturation point shifts to lower density and has a smaller 
binding energy per nucleon. In analogy with their result, we expect that a proper Dirac- 
Brueckner calculation for strange nuclear matter will also show a shift of saturation points 
as compared to what we find in our present Brueckner-Hartree-Fock calculations. However, 
we believe that the calculations presented here suffice for our purpose, which is to study the 
general features of strange nuclear matter. Any shift of saturation points is only expected 
to be relevant in those cases where the matter under consideration is on the boundary of 
being bound or unbound. In those cases, the Dirac-Brueckner result might show that matter 
which we find to be just bound is actually unbound, or vice versa. 

A second approximation is that the single-particle potential is radically put to zero for 
momenta Pi above the Fermi sea: the so-called "standard" choice. This causes a discontin- 
uous jump in Eiipi) at Pi = pp, and so it is also known as the "gap" choice. An alternative 
choice is to retain a nonzero value for momenta above the Fermi sea: the so-called "contin- 



uous" choice [23]. There are various physical arguments which favor this latter choice [23 



but its main effect is to merely shift the saturation curve to give more binding, without 



changing the overall density dependence very much; see, e.g., Refs. [0,0 for the effect in 
ordinary nuclear matter. However, these differences are only of relevance on a quantitative 
level (e.g., when a comparison is made with the experimental saturation point), and so we 
argue that for this first study of the general features of strange nuclear matter it suffices to 
use the gap choice. Another motivation is that the continuous choice considerably compli- 
cates the propagator in the Bethe-Goldstone equation, which makes the calculations much 
more cumbersome and computer intensive. 

A further simplification in the self-consistency calculation is that we assume a quadratic 
momentum dependence for the single-particle potential U(p). This means that we can define 
an effective baryon mass M* in terms of which the single-particle energy can be written as 

2 

E l (p) = M l + ^ + neU l (p) (5a) 
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2 

M l + J^ + Ileum, (5b) 



where Ui(0) is easily obtained from Eq. (|]), while M* is obtained from 

(6) 



Ml 



l+ TZeU i ipp)-TleU i (0) 1 1 



p 2 F /2M t 



The advantage of using Eq. ( |5"fj| ) is that the self-consistency condition only needs to be cal- 
culated at p = and p = pp, rather than at a range of momentum values < p < pp, as 
required when using Eq. (|5a|). Also, the binding energy does not require a numerical integra- 
tion, but is easily done analytically. We checked for various cases that the parameterization 
of Eq. ( |5bD indeed fairly accurately represents the single-particle energy as obtained from 



an explicit calculation using Eq. (5"aQ 



In Eq. (|2|), the Pauli operator needs to be expressed in terms of P and q. Clearly, 
|jLi]P+q| > kp and I^P— q| > kp are both satisfied when q > fiiP+kft and q > ^P+kp . 
Similarly, when q 2 < k^p 2 — (j^iP) 2 or q 2 < kp^ 2 — (/i 2 P) 2 at least one of the inequalities is not 
satisfied. For values of q between these two limits, there are two restrictions on the angle 
6>(P, q), namely 

cos > — cos 1 = , 

2\x x Pq 

^ a- (/^) 2 + ? 2 ~ k F 2 m 
costf<costf 2 = — . (7) 

Since the angle #(P,q) is integrated over, we can approximate this latter constraint by 
taking an average value for the Pauli operator Q. We therefore define 

Q{q,P) = l, if q > max[ Ml P + k? , fi 2 P + k {2) ] 

= 0, if q 2 < max[4 1)2 -( Ml P) 2 ,4 2)2 -(// 2 P) 2 ] 

= min [cos 6*i, cos 6*2], otherwise. (8) 

In the partial-wave projection, the Bethe-Goldstone equation for a system with isospin 
T and total angular momentum J becomes 

" „"2 



G J Z )ls {q', q; P, 00) = Vff^tf, q) + ~ £ fdq" q 

7r frr„ J 



l"s h 

xvfipM, q") — r? ™ P l, n — ^ — ~ G i^ is(q", ?; p, <*>), (9) 

is, i s w ,h ' u0 _M-P^/2M-q" 2 /2^-X + ie < s > ZsW '^' ' hyi 

where X = for the gap choice and X = U(\fjLiP+q\)+U(\]j,zP— q|) plus angle- averaging for 
the continuous choice. The single-particle potential is obtained self-consistently from 

u(pi) = E t TW t W 271 r dcosd [** ^** G %.l k > h > p > E M + frfoOU 

2V#,s + 1 )( 2 *1 + X ) 1/0 ( 27F ) 

(10) 
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where we have explicitly separated off the angle dependence of the d 3 p 2 integral. The isospin 
factors are present to account for all the contributions of the possible isospin states. (Our 
calculations are done on the isospin basis.) Finally, the subscript AS denotes that we have to 
include both direct and exchanged (Hartree and Fock) contributions. For identical particles, 
this can be accounted for by multiplying the G matrix from the Bethe-Goldstone equation 
with the factor 1 — (—l) l+s+t , with t — 1 for singlet-even and triplet-odd partial waves and 
t = for singlet-odd and triplet-even partial waves. (Note that i is equivalent to the isospin 
in the case of pure NN or pure SS systems.) 

If we now want to include all species of the octet baryons, the above expressions can 
be easily generalized. First, the internal sum over I", s" in the Bethe-Goldstone equation 
then also involves a sum over all possible channels allowed for a particular two-baryon 
interaction. Of course, the propagator needs to be modified to account for the relevant 
masses and thresholds in each particular channel, and the Pauli operator should contain 
the Fermi momenta belonging to the relevant species. Second, the single-particle potentials 
have to be summed over all baryon species. Using the notation U B for the single-particle 
potential of particle B due to the interactions with particles B' in the medium, and a bra-ket 
notation for the final-initial state particles, we have 



pi), 



(ii) 



where 



_ (2J+l)(2T+i; 
tm,s (2«B + l)(2f B 



2tt / dcos6 



p^dp 2 



1) 7_i J (2tt 

x An(BB'\Gi T ls [k,k;P,E[ B \ Pl 



Ef\ V2 )} A s\BB'). 



(12) 



The allowed values of T,J,l,s depend on what particular baryons make up the scattering 
process B + B' — ► B + B' . For example, C/^ gets contributions from direct isospin-3/2 
EiV scattering, but also from the coupled-channel isospin-1/2 (AN, UN) scattering. In our 
calculations we include all partial waves up to J = 4. Finally, the binding energy per baryon 
is obtained from 



A B 



1 



-l)(2t B + l 
KeU B (0) 



(4 S) ) d 3 k 



(27T) 




(B)2 



(13) 



IV. RESULTS 



A. Pure systems 



The first interesting question to ask is whether all of the isospin symmetric matter 
(T z = 0) made of only one kind of hadrons is bound. Our results are displayed in Fig. for 
models NSC97e (dashed curves) and NSC97f (solid curves). 
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The saturation curves for the purely nuclear system are very similar to what is obtained 
for other one-boson-exchange AW potentials found in the literature. The results for the two 
NSC97 models are practically indistinguishable, which, in fact, is true for all six NSC97 
models. This is a reflection of the fact that these models all describe the AW scattering data 
equally well. 

We see that the pure A system is not bound at all for both models. The pure A system was 
calculated including the coupling to the EN and EE channels. If this coupling is switched 
off (and so only elastic AA scattering is possible), the curves shift to slightly higher values: 
about 10% higher for NSC97e and about 5% higher for NSC97f. 

The pure H system is more bound and saturates at higher densities than the nucleon 
system. The difference in the results for NSC97e and NSC97f is due to the fact that the 
3 Si partial-wave contribution to the single-particle potential (which is large and positive) for 
NSC97f is almost 40% larger than for NSC97e. In addition, the 3 P2 partial- wave contribution 
(which is also large, but negative) for NSC97f is about 10% less attractive and largely 
compensates for an increased attraction in the 1 Sq partial- wave contribution. An important 
part of the attraction is due to the scalar-exchange part of the potential. Since the existence 
of a nonet of scalar mesons with masses below 1 GeV/c is still highly controversial (especially 
the low-mass isoscalar a meson), a comment is in order. We first note that within a one- 
boson-exchange model for the AW interaction, the scalar-exchange contribution plays a 
crucial role in providing the required attraction. Whether this contribution represents a 
true exchange of scalar mesons or just an effective parameterization of two-pion exchange 
and more complicated interactions is a question which goes beyond the scope of this paper. 
Here we only want to mention that arbitrarily removing the a (or, in our case |I9 |, the broad 



e) contribution renders purely nuclear matter unbound at all densities. However, the pure 
S system still remains bound, although it saturates at a smaller density. 

Pure £ matter with T z = is expected to be highly unstable, because it can decay 
strongly into A matter via — > AA and — > AA. However, we can consider pure 
S + or pure E~ matter with T z = ±2, respectively. Inclusion of the Coulomb interaction in 
this charged system is beyond our present calculation (and is likely to modify the result), 
and so the third panel of Fig. |2] represents the result without the Coulomb interaction and is 
included for illustrative purposes only. The difference in the results for NSC97e and NSC97f 
is due to the fact that the single-particle potential for NSC97f is slightly more attractive for 
almost all partial-wave contributions, which adds up to a substantial difference in the total 
single-particle potential. 

The results of these calculations show that for the employed baryon-baryon models of 
Ref. ]r9"| , |2"0"! the only possible long-lived strange pure system within SU(3) is the H system, 
which is stable against strong decays. This suggests that the inner core of high-density 
astrophysical objects could be rich in H particles. 



B. {N, A} systems 

We next consider the change of the nuclear binding by adding A's. The results are 
shown in Fig. |3|. We should point out that, even though the £ and H are not explicitly 
included as part of the medium, the coupling to these particles via the transition potentials 
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V(AN -> EAT) and ^ (AA -> HA", EE) are included. We see that the binding first slightly 
increases until the A fraction reaches about 10%. The system rapidly becomes less bound 
when more A particles are added. 

The most stable system occurs at about p = 0.28 fm~ 3 and xa — 0.1 with E/Ab = —13.2 
MeV for NSC97e and at about p = 0.27 far 3 and xa = 0.05 with E/A B = -12.7 MeV 
for NSC97f. The A fraction at the minimum is very similar to the Xa ~ 0.06 found in 
Ref. [ [L2[| , but in their case the minimum occurs at a lower density with a higher binding 
energy (p = 0.23 fm~ 3 with E/Ab — —16.0 MeV). This difference is due to various reasons. 
First of all, the larger binding energy found in Ref. [EJ is due to the fact that they use 
the continuous choice for the single-particle energy, whereas we use the gap choice. The 
smaller density at which their saturation occurs is due to the different AW potential that 
is employed: they use the parameterized Paris AW potential f2"B| , whereas here we use the 



AW potential as given by the NSC97 models. Furthermore, the results of Ref. |12| do not 
include the effect of A's interacting with themselves. We do include the AA interaction 
and, although we find that the contribution of U^' to the total single-particle potential is 
rather small, its influence starts to become noticeable as the A density increases. At low A 
density the individual partial-wave contributions to are almost all negative, but as the 
A density increases some of them start to give (relatively important) positive contributions. 
As a consequence, the effect of the inclusion of the AA interaction for increasing xa is to 
shift the density at which saturation occurs to lower values. The effect is most pronounced 
for NSC97f. In Ref. Jl2 the saturation density is more or less independent of the A fraction; 
see their Fig. 5. 

Our results suggest that multi-A systems produced in relativistic heavy-ion collisions 
through, for example, the coalescent mechanism could be loosely bound. However, the 
presence of a large fraction of A's in the inner core of a dense star seems unlikely, since a 
too large fraction tends to destabilize it. 



C. {TV, A, S} systems 



A further difference between the present work and that of Ref. |T2|] is that here we can also 
include the 5 (and E) as part of the medium. To investigate the influence of including H's, 
we perform calculations for a system consisting of N, A, and H. The E's are excluded since 
they can easily be annihilated, as stated before. Another motivation for their exclusion 
is that the Q-values for the strong transitions EA" — > AA^, EE — ► AA, EA — > EN, and 
EH — > AH are about 78, 156, 50, and 80 MeV, respectively. To Pauli block these processes, 
we need a rather high density of A. On the other hand, the Q-value of HA" — > AA is only 
about 28 MeV. The presence of H could then help prevent the collapse of the {A^, A,H} 
system since the AA — > HA" reaction can be Pauli blocked. The importance of including H's 
was first pointed out in Ref. ||. 

It is interesting to first investigate the charge-neutral systems consisting of only N, A, 
and H. They can be formed with a density distribution of p p = p n = p= = Ps~- Our results 
are shown in Fig. |j. We see that the systems are only loosely bound. The A density is too low 
to Pauli block the HA" — > AA process. This suggests that a charge-neutral strangeness-rich 
system is unlikely to be seen in nature or to be created in relativistic heavy-ion collisions. 
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We now turn to investigating the dependence of the binding energy of the N, A, S system 
on the strangeness per baryon. As discussed in previous works @-ll], this dependence is 



most relevant to the investigation of relativistic heavy-ion collisions. For this purpose, it is 
useful to define the fractions x% — Pi/ P f° r the different species i = N,A,E, as in Eq. (|J). 
The strangeness per baryon for the {N, A,H} systems can then easily be calculated as 
fs — Xa + 2%=. Since we work on the isospin basis, we will only consider systems with 
Pn = P P = \pn and p=o = p s - = \pu- The charge per baryon is then f q = |(xjv — Xs)- For 
each {N, A, 5} system with a given xn we carried out calculations for various combinations 
of (xa + Xs) = (1 _ Xn)- I n general, we find that the system is less bound for high Xs, as 
could already be inferred from comparing Figs. and [|. The reason is that ui E) becomes 
more negative (attractive) as xs increases, but this is compensated by positive (repulsive) 
contributions from Upf A and U S N . The cancellations are large enough to prevent p^ from 
becoming too large, and so the large binding energies as found close to the saturation point 
of the pure S system (see Fig. ||) cannot be reached. 

These calculations allow us to examine the f s (f q ) dependence of the binding energies 
at the saturation point of an {A^, A, S} system. The results are shown in Fig. [|. In each 
case, the curves cover the allowed f s values that can be reached. We see that as f s (and 
f q ) increases, the system becomes less bound. When xn = 0.7, the purely {A^, A} system 
(i.e., xs — 0) has the lowest binding, but when xn gets smaller than 0.6, systems with 
increasing xs are preferred. This follows from the fact that in those cases the curves show 
a (shallow) minimum. However, when xn ~ 0.4 the system becomes unbound. Our results 
contradict the results from the relativistic mean-field calculations of Refs. PHTTf (see Figs. 2 
and 3 of Ref. |§). The differences, of course, could be due to finite-size effects, since these 
authors consider a shell model for finite {A^, A, S} systems (up to very large = 310). 
But it is more likely that the difference is due to the differences in dynamical content of the 
calculations. This can be understood by observing that the mean-field calculation within 
the Brueckner-Hartree-Fock approach amounts to neglecting the residual baryon-baryon 
interaction terms in calculating E/Ab- Namely, the mean- field results can be obtained from 
Eq. ([I"3| ) by making the change \Us{k) — > Usik) in Eq. (|T3"D. In all of the cases, Usik) is 
comparable to kinetic energies and, hence, the finite-size mean-field results can be radically 
different from our full calculations. 



D. Inclusion of 5] 

Although the S particles can be easily annihilated, as stated before, their presence largely 
increases the binding of the system, which might be of relevance in the formation and stability 
of such a system. To demonstrate the effect of adding S's we first consider a charge- neutral 
{Af, A,S} system with xn = Xs, and so f s = 1. For each choice of xn = Xh we can add 
different fractions of S's as long as (xa + Xs) = (1 — Xn — Xs)- Note that these systems 
all still have f s = 1. In order to prevent having to show numerous figures or tables, we 
will here only consider the systems with xa — Xs- The results are shown in Fig. |6]. They 
should be compared with the saturation minima in Fig. |J We clearly see that adding also 
S's drastically improves the binding of the system compared to adding only A's. 

The increase in binding due to adding S's allows us to further explore the f s dependence 
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in the region beyond f s < 1. Again, numerous choices for the different particle fractions 
are possible, but here we will only restrict ourselves to systems with equal fractions for the 
strange particles, i.e., xy = Xa = Xs — Xs- We find that these systems are bound for any 
value of xy- The f s dependence of this mixed system is shown in Fig. [7[ We see that as f s 
increases, the system becomes more bound, except in the low f s region. Note that at the 
highest f s = 2, the system is the pure H system as in Fig. |[ Of course, our special choice for 
Xy means that the strangeness per baryon can only go up to f s = 4/3, and so the result for 
f s > 4/3 in Fig. [7] is not calculated, but is rather obtained by simply extrapolating to the pure 
H result which has / s = 2. Comparing with the results of Fig. ^ for the {N, A, 5} system, 
it is clear that the addition of the S component drastically changes the f s dependence, but 
it is important to remember here that the system containing £ is highly unstable against 
strong decays, as discussed above. The results shown in Fig. [?| perhaps cannot be verified 
in relativistic heavy-ion collisions. However, it could represent the situation in the core of 
neutron stars in which the presence of high density e~ can produce a lot of £~, initiated by 
the reaction e~ + n — > + v. 



V. CONCLUSION 

We have investigated strange hadron matter in the context of a baryon-baryon potential 
model based on SU(3) symmetry. The parameters of the potential model were fitted to 
the NN and YN scattering data, and the assumption of (broken) SU(3) symmetry allows 
us to extend the model to also describe the other interaction channels that are allowed 
for the baryon octet; i.e., the YY, EN, HY, and SS interactions. The potential for these 
interactions is defined without the necessity of having to introduce new free parameters. 

The calculations have been carried out by using the Brueckner-Hartree-Fock approxima- 
tion. Within the framework that we use to define the potential model for the baryon-baryon 
interactions, we find that the pure A system is unbound, whereas the pure S system is more 
strongly bound and saturates at higher densities than the pure N system. Adding A's to 
pure nuclear matter slightly increases binding, as long as the A fraction is less than about 
10%. Larger fractions cause a decrease in binding. Adding S's, of importance due to the 
reaction AA — > EN, drastically reduces binding, and so {N, A, E} systems are only (weakly) 
bound for nucleon fractions larger than 40%. Our results represent a step forward, since the 



previous Brueckner-Hartree-Fock calculation 12] did not include the S's 



By carrying out extensive calculations for the {N, A, E} and {A^, A, £,5} systems, we 
have predicted the dependence of the binding energy on the strangeness per baryon, a quan- 
tity that is needed to be determined as precisely as possible for identifying the strange quark 
matter created in relativistic heavy-ion collisions. Our results are significantly different from 
previous calculations based on relativistic mean-field models. We argue that the differences 
are mainly due to the two-body correlations, which are neglected in relativistic mean-field 
models. 

To close, we would like to point out that the present work is merely a first step towards 
a rigorous many-body calculation. In the future, we need to investigate the three-body 
terms in the Brueckner-Hartree-Fock approach. This is particularly important if one is to 
give more precise predictions in the higher-density regions. Another reason is that perhaps 
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it is only at that order that the issue of using the gap choice or the continuous choice for 
the single-particle energies in solving the self-consistent Eqs. (§) and (U) becomes irrelevant, 
as was shown in a recent work on ordinary nuclear matter P7j. Our investigation in this 
direction will be published elsewhere. 

ACKNOWLEDGMENTS 

This work was partly supported by the U.S. Department of Energy, Nuclear Physics 
Division, under Contract No. W-31-109-ENG-38. 



12 



REFERENCES 



[1] As reviewed by B. Shiva Kumar, Nucl. Phys. A590, 29c (1995). 

[2] As reviewed by G. Baym, Nucl. Phys. A590, 233c (1995); M. Prakash and J.M. Lat- 

timer, Nucl. Phys. A639, 433c (1998), and references therein. 
[3] A.R. Bodmer, Phys. Rev. D 4, 1601 (1971). 

[4] S.A. Chin and A.K. Kerman, Phys. Rev. Lett. 43, 1292 (1979); E.P Gilson and R.L. 

Jaffe, Phys. Rev. Lett. 71, 332 (1993). 
[5] E. Witten, Phys. Rev. D 30, 272 (1984). 

[6] E. Farhi and R.L. Jaffe, Phys. Rev. D 30, 2379 (1984); 32, 2452 (1985). 

[7] J. Schaffner-Bielich, C. Greiner, A. Diener, and H. Stocker, Phys. Rev. C 55, 3038 

(1997) . 

[8] V.R. Pandharipande, Nucl. Phys. A178, 123 (1971). 

[9] J. Schaffner, C.B. Dover, A. Gal, C. Greiner, and H. Stocker, Phys. Rev. Lett. 71, 1328 
(1993). 

[10] J. Schaffner, C.B. Dover, A. Gal, D.J. Millener, C. Greiner, and H. Stocker, Ann. Phys. 

(N.Y.) 235, 35 (1994). 
[11] A. Gal and C.B. Dover, Nucl. Phys. A585, lc (1995). 

[12] H.-J. Schulze, M. Baldo, U. Lombardo, J. Cugnon, and A. Lejeune, Phys. Rev. C 57, 
704 (1998). 

[13] C. Greiner, P. Koch, and H. StScker, Phys. Rev. Lett. 58, 1825 (1987). 

[14] C. Greiner, D.-H. Rischke, H. StScker, and P. Koch, Phys. Rev. D 38, 2797 (1988). 

[15] C. Greiner and H. St6cker, Phys. Rev. D 44, 3517 (1991). 

[16] H.J. Crawford, M.S. Desai, and G.L. Shaw, Phys. Rev. D 45, 857 (1992). 

[17] A.J. Baltz, C.B. Dover, S.H. Kahana, Y. Pang, T.J. Schlagel, and E. Schnedermann, 

Phys. Lett. B 325, 7 (1994). 
[18] P.M.M. Maessen, Th.A. Rijken, and J.J. de Swart, Phys. Rev. C 40, 2226 (1989). 
[19] Th.A. Rijken, V.G.J. Stoks, and Y. Yamamoto, Phys. Rev. C 59, 21 (1999). 
[20] V.G.J. Stoks and Th.A. Rijken, submitted to Phys. Rev. C. 
[21] K.A. Brueckner and J.L. Gammel, Phys. Rev. 109, 1023 (1958). 
[22] R. Brockmann and R. Machleidt, Phys. Rev. C 42, 1965 (1990). 
[23] J. -P. Jeukenne, A. Lejeune, and C. Mahaux, Phys. Rep. 25C, 83 (1976). 
[24] M. Baldo, I. Bombaci, L.S. Ferreira, G. Giansiracusa, and U. Lombardo, Phys. Rev. C 

43, 2605 (1991); H.Q. Song, M. Baldo, G. Giansiracusa, and U. Lombardo, Phys. Lett. 

B 411, 237 (1997). 

[25] H.-J. Schulze, J. Cugnon, A. Lejeune, M. Baldo, and U. Lombardo, Phys. Rev. C 52, 
2785 (1995). 

[26] M. Lacombe, B. Loiseau, J.M. Richard, R. Vinh Mau, J. Cote, P. Pires, and R. de 

Tourreil, Phys. Rev. C 21, 861 (1980). 
[27] H.Q. Song, M. Baldo, G. Giansiracusa, and U. Lombardo, Phys. Rev. Lett. 81, 1584 

(1998) . 



13 



FIGURES 

FIG. 1. 1 S'o phase shifts for elastic identical-particle scattering, for models NSC97e and NSC97f. 

FIG. 2. Saturation of pure systems. The dashed curve and solid curve are the predictions 
of model (e) and (f), respectively. The pure £ system represents the T z = ±2 case without the 
Coulomb interaction. 

FIG. 3. Saturation of {N, A} systems for various fractions pa/(pn + Pa)- 

FIG. 4. Saturation of charge-neutral {N, A, E} systems for various fractions pa/(pn + Pa + ps)- 

FIG. 5. f s dependence of the {A^, A, E} system for three given fractions of nucleons xn- 

FIG. 6. Saturation of charge neutral {A^, A, system with \N = Xs an d XA = Xs, as a 

function of XA+S- 

FIG. 7. f s dependence of the {N, A, S, E} system with xa = XT, = Xs- 
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